Force balance and membrane shedding at the Red Blood Cell surface 
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During the aging of the red-blood cell, or under conditions of extreme echinocytosis, membrane 
is shed from the cell plasma membrane in the form of nano- vesicles. We propose that this process is 
the result of the self- adapt at ion of the membrane surface area to the elastic stress imposed by the 
spectrin cytoskeleton, via the local buckling of membrane under increasing cytoskeleton stiffness. 
This model introduces the concept of force balance as a regulatory process at the cell membrane, 
and quantitatively reproduces the rate of area loss in aging red-blood cells. 



The membrane of red-blood cells (RBCs) is a fluid 
sheet of lipids and proteins, connected through node com- 
plexes to a two-dimensional cytoskeleton primarily com- 
posed of flexible spectrin protein filaments Q . The me- 
chanical properties and the shape of the RBC result from 
a competition between the bending elasticity of the cell 
membrane and the shear and dilation elasticity of the 
cytoskeleton 0, Because of the relative structural 
simplicity of the spectrin network, the RBC is a system 
of choice to study the co-organization of the lipid bilayer 
and the cytoskeleton cortex at the plasma membrane [2j • 

Over the natural life-span of the RBC, the cell looses 
a large fraction of its membrane (~ 20% over 120 days 
in humans and ~ 10% over 50 days in rabbits 0,0) by 
the shedding of small vesicles containing hemoglobin but 
devoid of cytoskeleton 0, IS • The shed vesicles have a 
size (100-200nm Q) similar to the mesh-size of the un- 
derlying spectrin network an d mav take the form 
of elongated cylinders . The shedding of vesicles is also 
observed under conditions that trigger the transition into 
echinocyte shape, such as ATP depletion and Ca 2+ load- 
ing In this case, the shedding of cytoskeleton-free 
vesicles often occurs at the tip of membrane spicules con- 
taining cytoskeleton 0, 0, [3] • The shedding of mem- 
brane area is crucial to the fate of the cell, as a decreasing 
surface-to- volume ratio is thought to be connected to the 
phagocytosis of old RBCs. It is this process of membrane 
shedding that we model in this work. 

When the entire membrane of the RBC is removed un- 
der physiological conditions, the cytoskeleton shrinks to 
a 3 to 5-fold smaller area [13j ■ This indicates that the cy- 
toskeleton is stretched by its attachment to the cell mem- 
brane 0, 0] , and therefore exerts a compression force on 
the membrane. As the cell ages, its observed stiffness 
increases by -20% [lEIH, 

as does the density of its cy- 
toskeleton network by about — 30 — 40% 17]. A recent 
model 18] proposes that these changes may be related to 
the diminishing levels of ATP as the cell ages jjjij) . This 
model relates the ATP level to the degree of dissocia- 



tions in the cytoskeleton network; Lower ATP levels cor- 
respond to a stronger cytoskeleton network, and results 
in larger compressive forces on the cell membrane. These 
forces are balanced by the membrane bending stress [jj, 
which may either be of entropic origin (from constraints 
on the membrane thermal fluctuations), or resulting from 
an overall membrane curvature (blobbing). We show be- 
low that beyond a compression threshold, the membrane 
undergoes a buckling transition at the scale of the spec- 
trin mesh-size, and one or a few membrane blebs grow 
in size in order to release the overall stress. We first 
analyze the buckling transition of a single membrane- 
cytoskeleton unit. The thermally fluctuating membrane 
is assumed flat at large lengthscales (we do not account 
for the global shape of the RBC), while the spectrin fila- 
ments are treated as phantom Hookean springs |20|| . at- 
tached at their ends to the membrane (FigOJ. We then 
consider a large elastic network attached to a membrane, 
and show that the lowest energy state beyond a critical 
compression strain consists of an essentially stress-free 
membrane connected to a large bleb which incorporates 
all the excess membrane area. We argue that such a state 
is never reached because the bleb pinches off the mem- 
brane to form a free vesicle before reaching large sizes. 
Finally, we relate our results to the observed shedding of 
vesicles during RBC aging. 



Let's assume that one cytoskeleton unit, of area a, is 
connected to a patch of membrane of area Strue (Fig^l, 
and constrains its fluctuations |2l| . There exist an opti- 
mal unit area that leaves the fluctuations unconstrained. 
This area, which we call s, is the apparent area of a freely 
fluctuating membrane and is slightly smaller than the 
true area by a factor s/s tr ue = (1 — l/adogs trMe /Z 2 ) [22]] . 
Here, a = 8irn/k B T is the bending energy of a closed 
sphere (compared to the thermal energy), and I a small 
size of order the membrane thickness. If the membrane 
is stretched (a > s), or mildly compressed, it remains flat 
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FIG. 1: Energy of a patch of fluctuating membrane of ap- 
parent area s fit to an elastic frame of rest area ao. The 
energy shows an inflexion point when the membrane buck- 
les. The gray line is the numerical optimization of Eq.© 
with Eq.(|5J and the dashed line is the asymptotes Eq.JlJ, for 
a = 500 (k = 20fc B T) and e = 5 10" 2 (e = 10" 5 J/m 2 and 
ao = (lOOnm) 2 ). Inset: The deformation of the cytoskeleton 
(area a) is maximal at the transition. 



on average and its free energy is |22|, |2 



ffl[s,a] 



-k R T 



(e Q (°/s-i) - 1 - a (a/s - 1)) (1) 



This energy is quadratic (~ Tr/2a 2 (a/s — l) 2 ) for small 
compression (a(s/a — 1) < 1), and becomes linear (~ 
7r/2a(l — a/s)) for large compression. The energy of the 
flat compressed patch eventually becomes larger than the 
bending energy of a spherical cap, leading to the breaking 
of the up/down symmetry and the formation of a curved 
bleb. We investigate the buckling transition by assuming 
that a buckled membrane fluctuates around a mean shape 
defined by a spherical cap of area s& > a, to be optimized 
(Fig.l). The energy of a buckled patch reads 



fm[s,a] 



Bttk 1 



ffl[s,s b ] 



(2) 



The first term is the bending energy of a spherical cap of 
area Sb on a frame of area a, and the second the energy 
of the thermal fluctuations (Eq.QJ), assumed unaffected 
by the curvature of the bleb. The membrane is flat on 
average if sj, = a and is buckled if Sb > a. Comparing 
Eq.Q and Eq.JSJ shows that in the relevant limit (a ^> 
1), the optimal buckled area is Sb — s — 2a/ira, and 
buckling is to be expected for s/a > 1 + 2/ira. 

In addition to the membrane entropic elasticity f m , 
the total energy of a membrane unit has a contribution 
from the cytoskeleton elasticity, treated here as a simple 



/ = 77 ea o ( — - 1 ) + fm [s, a] 
2 \a 



(3) 



Here, e and ao are the stretching modulus and rest area 
of the cytoskeleton unit mesh. The full free energy for 
one patch (Eq.©, with Eq.J2J)) is to be optimized with 
respect to two areas: the buckled membrane area s&, and 
the area of the stretched cytoskeleton mesh (a > ao). The 
buckling transition is controlled by the area mismatch 
s/ao and the ratio of the typical cytoskeleton to mem- 
brane elastic parameters e = eao/ak B T. In the case of 
the RBC, we expect e«l and the buckling transition is 
second order (the cytoskeleton area is continuous through 
the transition, Fig^- inset). The numerical minimization 
of the full energy in shown in Fig^ The buckling tran- 
sition occurs for a critical mismatch s/ao — l/Ve (for 
a»l and e«l), and the optimal meshsize and energy 
below and above the transition (subscripts "/" and "b" , 
respectively) are: 
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Note that the area of the stretched cytoskeleton mesh in- 
creases with increasing membrane area s below the buck- 
ling transition, but decreases with increasing s above the 
transition (Inset of Fig^l. 

We now proceed to describe the buckling instability 
of many connected patches. Under low compression, the 
membrane of the cell is equally shared by all cytoskele- 
ton units, and if the membrane were not fluid, an increase 
of the cytoskeleton compressive stress would eventually 
lead to the buckling of every single membrane patch. Bi- 
ological membranes are fluid, so once a unit buckles, it 
may grow into a large bleb by feeding from the excess 
area of the surrounding units. This reduces their energy 
without much increase of the energy of the bud, which 
saturates to the energy of a full sphere a for large area 
(Eq.Q). Mathematically, coexistence between flat and 
buckled units is allowed by the existence of an inflex- 
ion point in the energy of one isolated membrane patch 
(FigD). 

Let's consider N membrane-cytoskeleton units similar 
to the one shown in Fig^ with total membrane area S. 
We investigate a membrane composed of (N — 1) units 
of area s that are in the flat (unbuckled) state, plus one 
unit of area S — (N — l)s, that may be buckled if such 
shape lowers the total energy: 

Ffl = (N- l)f f (s) + f f . b (S-(N- 1) 8 ) (5) 

where the energies in the flat and buckled state (// and 
fb respectively) are given in Eq.Q. The parameter that 
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FIG. 2: Energy of a collection of iV membrane + cytoskele- 
ton units as a function of the membrane area per unit s (from 
Eq.@). Darker shades of gray corresponds to larger total 
membrane area S. The energy shows a narrow local mini- 
mum for an homogeneous flat membrane at s — S/N. A local 
minimum at s ~ ao appears beyond the binodal point. It cor- 
responds to stress-free units coexisting with a large membrane 
bleb. The homogeneous membrane (a) is globally unstable for 
S > S\ ■ After a waiting time depending on the nearness to 
the spinodal line (Eq.JJJ) a small bleb forms (b), grows (c), 
and stochastically pinches off the cell (d) reducing the total 
membrane area and leaving the system in a deeper local min- 



imum until the next blebbing event. N 
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and S/Nao — 1.01 (binodal), 1.063 (equilibrium), 1.085, and 
1.095. The spinodal point {S/Na = 4.98) is not shown. 



controls the state of the membrane is the dimensionless 
average area per unit (area mismatch) S/Nao, and the 
state of the membrane is unambiguously given by the 
area of the flat units s. For small area mismatch, the 
only equilibrium state corresponds to a uniform flat mem- 
brane with s = S/N. Upon increase of the mismatch, 
the system eventually reaches the binodal point. A local 
minimum of the total energy F tot appears for s < S/N, 
corresponding to an inhomogeneous membrane state with 
one buckled unit (FigJ2J. For large area mismatch, the 
optimal conformation is such that all the flat units are 
near their optimal area (s > ao), with all the excess area 
concentrated in the one nearly spherical bleb of energy 
F b tot ~ a. 

The coexistence between flat and buckled states re- 
quires the equality of tensions: d s ff(s) = d s fb(S — (N — 
l)s) (-> d s F b tot = 0). The buckled state is the equi- 
librium state beyond a total area S\ for which buck- 
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led and flat states have the same energy F b 
Nae/2(S\ aq /(Na ) - l) 2 . Although it is thermodynami- 
cally favorable to form a large bleb devoid of cytoskeleton 
for S > S\ , this process requires an activation energy 
AF (see Fig|5J). The energy barrier is very large (of order 
a ~ 500fc B T) close to equilibrium, decreases with increas- 
ing area mismatch, and vanishes at the "spinodal" point 
Si in for which every single unit is unstable with respect 
to buckling. In the relevant limit (e <C 1 and iVe ^> 1), 
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the equilibrium and spinodal mismatches are given by: 

S | 1 

' Na lspm -^ 

After the buckling of one membrane unit (Fig|2b), the 
system then evolves downhill (Fig|2t) toward its global 
equilibrium state, which corresponds to a large bleb in 
coexistence with stress-free units (s ~ ao). The "equilib- 
rium bleb" takes all the excess area, and should thus be 
very large (area of order aoy/N/e if buckling occurs at 
equilibrium and agN/y/E near the spinodal curve). 

Although such large blebs are observed in blebbing 
cells, they generally involve the unbinding of large 
patches (~ /^m 2 ) of membrane from the cytoskeleton |24| 
and don't bud off the cell. In aging RBC, small spherical 
or cylindrical vesicles of size comparable to the cytoskele- 
ton meshsize (~ (lOOnrn) 2 ) are being shed. This phe- 
nomenon could have a kinetic origin and be explained by 
the following simple argument. Once a unit membrane 
has buckled, it initially grows to form a hemispherical 
cap (of bending energy a/2), from which it may grow 
further into a shape approaching either a fully formed 
sphere or a cylinder. The bending energy of the cylinder 
(~ a(s&/(16a) + 3/8) is smaller than the energy of the 
sphere (~ a(l — a/sb)), as long as s& < 8a, so we argue 
that the buckled unit should initially grow as a cylinder. 

The growing bleb should undergo large fluctuations, 
owing to the fact that the membrane is under small (or 
even negative) surface tension. One may thus expect the 
high curvature neck connecting the bleb to the cell mem- 
brane (FigEt) to eventually close and the bleb to pinch- 
off before reaching a large size, in agreement with the 
experimental observations. Vesicle scission is expected 
to be much faster than the RBC aging process and could 
be promoted by phase separation of lipids and proteins 
within the membrane of the bleb [23 ■ Once shedding has 
occurred, the system finds itself in a new (deeper) local 
energy minimum (Fig|21i), corresponding to the reduced 
total membrane area S. 

In the context of the RBC, the cytoskeleton stiffness e 
continuously increases with time due to the depletion of 
the ATP inside the RBC 18). Correspondingly, the spin- 
odal and equilibrium areas (Eq.JHJl) decrease with time, 
and so does the energy barrier that keeps the cell in a 
stressed state. The buckling rate should eventually ad- 
just to the rate of cytoskeleton stiffening, as the pinching- 
off of the buckled membrane reduces the RBC area and 
stabilizes the flat state until further stiffening. Blebbing 
occurs when a thermal fluctuation moves the system over 
the energy barrier. This requires a waiting time of order 
Toe AF / kBT , where tq ~ 0.1ms is the typical timescale 
for membrane fluctuations on the length-scales of a sin- 
gle cytoskeleton unit [23 ) an d AF is the area-dependent 
energy barrier (FigEJ. Accordingly, the RBC area at a 
given time is such that AF = fc B Tlogr e /ro, where r e is 
the characteristic time of cytoskeleton stiffening. 
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attached elastic network. 
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FIG. 3: "Non-equilibrium Phase diagram" for RBC buckling. 
Cells with large membrane area (y-axis) are unstable to buck- 
ling. The critical area for buckling depends upon the (time- 
dependent) cytoskeleton stiffness (x-axis), and the rate at 
which the cytoskeleton stiffens. Aging RBCs follow trajecto- 
ries showing cytoskeleton stiffening and the loss of membrane 
area (through the shedding of vesicles, represented as vertical 
arrows) . Trajectories are such that the rate of vesicle shedding 
matches the rate of cytoskeleton stiffening (r e = roe AF ' s ' <! ' 
see text). Three rates are represented (r e = lday, 1 second, 
and the spinodal line, Eq.JfjJ, corresponds to infinitely fast 
stiffening). Parameters are the same as Figs. 1-2 



Fig[3]shows the evolution of the cell membrane area for 
various rates of cytoskeleton stiffening. Buckling occurs 
close to S — S\ S pi n in the limit of a very fast ATP de- 
pletion rate, and at S = S\ eq for unrealistically slow rate 
(not shown). Near the spinodal line, an increase Se of the 
cytoskeleton stiffness leads to a decrease SS = —S8e/(2e) 
of membrane area (Eq.JJJJ). The observed rate of mem- 
brane shedding from an aging RBC is 20% in ~ 100 
days, which translates to about one lOOnm-sized vesi- 
cle an hour, consistent with an energy barrier at the 
transition AF ~ 15fc B T. We predict this loss of area to 
correlate with a 40% stiffening of the spectrin cytoskele- 
ton (FigEJl, in agreement with mechanical and structural 
measurements [lj, [la, ll7J • The calculated metastable ex- 
cess membrane area S/Nao ~ 4 (FiglSt also agrees well 
with experimental observation |l3| . 

We conclude that the compression forces produced by 
an elastic network on an attached fluid membrane may 
lead to buckling, and eventual vesiculation of the mem- 
brane. These conditions occur under slow RBC aging, 
during ATP depletion or Ca 2+ loading, which all result 
in cytoskeleton stiffening 0]. Note that Ca 2+ loading 
occurs regularly at the neural synapse, where internal 
vesicles respond by fusing and releasing their content [13] • 
Our work shows that a cortical network may act to me- 
chanically produces small vesicles, in a process entirely 
distinct from the classical Clathrin-based budding mech- 
anism, ft may therefore be relevant to some of the many 
biological systems containing fluid membranes with an 
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